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Tsallis has proposed a generalization of Boltzmann-Gibbs thermostatistics by introducing 
a family of generalized nonextensive entropy functionals with a single parameter q. These 
reduce to the extensive Boltzmann-Gibbs form for q = 1, but a remarkable number of 
statistical and thermodynamic properties have been shown to be (/-invariant - that is, valid 
for any q. In this paper, we address the question of whether or not the value of q for a 
given viscous, incompressible fluid can be ascertained solely by measurement of the fluid's 
hydrodynamic properties. We find that the hydrodynamic equations expressing conservation 
of mass and momentum are q- invariant, but that for conservation of energy is not. Moreover, 
we find that ratios of transport coefficients may also be g-dependent. These dependences 
may therefore be exploited to measure q experimentally. 
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I. INTRODUCTION 

A. Motivation and Historical Background 

The concept of extensivity is introduced early in most textbooks on thermodynamics and statistical physics. 
The requirement that the entropy be additive establishes the form of the Boltzmann-Gibbs distribution 
via a straightforward argument. Recently, Tsallis has proposed a generalization of Boltzmann-Gibbs 
thermostatistics by introducing a family of generalized entropy functionals with a single parameter q. The 
proposed generalization is best described by the following two axioms: 

Axiom 1 The entropy functional associated with a probability distribution f(z) is 



Axiom 2 The experimentally measured value of a phase function g(z) is given by the q-expectation value, 

G g [f]= [dz [f(2.)] q g(z). (2) 



From the first axiom, we note that S q [f] reduces to the Boltzmann-Gibbs entropy in the limit as q — > 1, 

Sx[f] = lim f dz {/(z) - [/(z)H = -k B f dz J(z)ln/(z), (3) 

a-* 1 q - l j j 

so the generalized thermostatistics includes the usual one as a special case. It differs most notably in that fact 
that neither the entropy S q [f] itself nor the observables G q [f] are extensive thermodynamic variables when 
q 1. In spite of this difference, a remarkable number of statistical and thermodynamic properties have been 
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shown to be g-invariant - that is, valid for any q whatsoever. These include the convexity of the entropy, the 
equiprobability of the microcanonical ensemble, the Legendre-transform structure of thermodynamics H, and 
Onsager Reciprocity ||, to name but a few. Other familiar properties, such as the Fluctuation-Dissipation 
Theorem, are not q- invariant, but have simple and straightforward generalizations to arbitrary q Q|. The 
implication is that the assumption of extensivity plays a role analogous to that of the parallel postulate of 
Euclidean geometry; one can deny it and still get perfectly self-consistent formulations of thermodynamics 
and statistical physics. 

Of course, all this would be but an idle (though undeniably interesting) mathematical exercise unless there 
were actual physical systems whose thermostatistics are best described by the generalized form with q 1. 
The exciting realization in recent years is that there do seem to be some of these. A very abbreviated list of 
examples is: 

• Stellar polytropes (e.g., globular clusters) were long known to possess kinetic equilibria for which there 
was no corresponding hydrodynamic variational principle until Plastino and Plastino Q showed that 
such a variational principle was possible only for q < 7/9. 

• Experimental studies of pure-electron plasmas in Penning traps have indicated that such plasmas 
turbulently relax to a radial density profile that does not maximize the Boltzmann-Gibbs entropy. It 
has recently been shown || that the observed profiles are consistent with q = 1/2, or perhaps slightly 
higher 0. 

• The ubiquity of Levy flights in physics can be explained P~^0[| by the fact that they are universal 
cumulative distributions - in the same sense that the Central Limit Theorem establishes the Gaussian 
as universal - arising from the generalized thermostatistics with q > 5/3. 

• Experimental observations of the velocity distribution of electrons undergoing inverse bremsstrahlung 
absorption give results consistent with q ^ 1 

Two natural questions arise at this point: What characteristics do physical systems with q 1 have in 
common? Is there a way to predict the value of q for a given physical system? Currently, there is more 
intuition about the first of these questions than the second. Systems that violate extensivity tend to do 
so because they have a long-range interaction potential, long-time memory effects, or a fractal space-time 
structure. The first of these qualities can make surface effects important, even in the thermodynamic limit. 
For example, it is straightforward to show that the total energy of a system of particles in D dimensions 
with interaction potential proportional to diverges if £ < D p^ j. The second quality can invalidate the 
Markovian assumption on which much of our physical intuition is based. The third quality can introduce 
scaling behavior with dimensionality not equal to that of the embedding space, and it can invalidate the 
Ergodic Hypothesis which also figures prominently in the justification of the Boltzmann-Gibbs distribution. 

B. Plan of this Paper 

While it is certainly comforting to find familiar properties of thermodynamics and statistical mechanics 
that are g-invariant - because this reinforces our existing intuition - it is no less important to clearly 
identify those features that are demonstrably not g-invariant. To abuse our above analogy with noneuclidean 
geometry, these are the phenomena that are the analog of the "triangular excess" of a polygon, or of the 
curvature tensor. The reason these features are interesting is that they are what allows us to experimentally 
distinguish systems with different values of q. At the level of kinetic theory, this is not difficult: The 
canonical ensemble distribution function is g-dependent, so its direct measurement p"H could provide a way 
to experimentally ascertain q. A more subtle question is whether or not different values of q can give rise to 
different hydrodynamic behavior. This is the central question that we address in this paper. 

We begin by developing and studying in some detail as simple a kinetic theory as we can imagine: We 
consider an ideal gas for general values of q. By "ideal" here, we mean only that the particles do not interact 
except in point collisions; we make no assumptions about the nature of those collisions. An objection that 
might be raised immediately is that an ideal gas is not the sort of system for which we would expect a 
violation of Boltzmann-Gibbs statistics, since it has no long-range potential. As mentioned above, however, 
there is currently no a priori way to determine q for a given physical system, so we are free to mandate that 
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q ^ 1 — at least as a mathematical exercise. We may suppose that there is some extenuating circumstance 
that somehow forces this ideal gas to have q ^ 1. For example, the particles might well carry a long-term 
memory of previous collisions, or their geometrical arrangement and/or the shape of their container might 
conspire to lead to gross violations of the Ergodic Hypothesis. In any case, there is some precedent for 
using this system to illustrate the application of the generalized thermostatistics: Plastino, Plastino and 
Tsallis |l4| considered the partition function and equilibrium properties of this very system, including the 
q-dependence of its specific heat. In this paper, we concentrate on the hydrodynamic behavior of this system. 

We next construct a kinetic theory for our general-q ideal gas. For the sake of simplicity, our kinetic theory 
is based on the Bhatnager-Gross-Krook (BGK) collision operator, which we generalize to arbitrary q, and 
for which we derive a g-invariant Jf-Theorem. Here it may be argued that the BGK collision operator is too 
naive to be used for our purposes, and we ought to have adopted a treatment based on the full Boltzmann 
equation. In defense of the BGK operator, however, we note that it is well known to produce the correct 
form of the viscous, compressible Navier-Stokes equations when q = 1, though the transport coefficients 
are different from those derived by the full Boltzmann treatment. The ratios of the transport coefficients, 
however, are more robust in this regard; for example, both the BGK and Boltzmann treatments for q = 1 
yield a ratio of bulk to shear viscosity of — 2/D, even though the absolute values of those viscosities are 
different. For these reasons, we stick to the BGK operator in this paper, focus our attention on robust 
results such as the form of the hydrodynamic equations and the ratios of the transport coefficients, and leave 
the more complicated Boltzmann analysis to future studies. 

We then derive the viscous, compressible hydrodynamic equations obeyed by the system, using a general- 
ization of Chapman and Enskog's asymptotic expansion in Knudsen number (the ratio of mean-free path to 
scale length) . These equations are the generalization to arbitrary q of the usual Navier-Stokes equations 
of hydrodynamics. We find that the Navier-Stokes equations expressing conservation of mass and momentum 
are q-invariant, but that for conservation of energy is not. Moreover, we find that ratios of transport coeffi- 
cients may also be g-dependent. These dependences may therefore be exploited to measure q by experiments 
at the level of hydrodynamics. Finally, in the process of our analysis, we show that q has a hard upper 
bound of 1 + 2/(D + 2) for systems of this sort. 



II. GENERALIZED HYDRODYNAMIC EQUILIBRIA 



A. Generalized Thermostatistics 



We first review the construction of the canonical ensemble distribution function using the generalized 
thermostatistics [Q. We maximize S q [f], given by Eq. (Q), subject to the preservation of various linear 
global functionals of /(z). By Tsallis' second postulate, Eq. (||), these are given by 

C\[f]= J dz [/(z)]V(z), (4) 

where the index i ranges from 1 to the number of conserved quantities n. We are thus led to the variational 
principle, 

= <j|s g [/]-£AiC$[/]j, (5) 

where the Aj's are Lagrange multipliers. It is an elementary exercise to verify that this yields the equilibrium 
distribution function 




The n constants \ are then determined by the n Eqs. (^) which may be written 
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7 '(z) 



In passing, we note that a very recently proposed modification to Tsallis' second axiom | 
the ^-expectation values as follows: 



(7) 



L 611 would normalize 



G'lf) 



Jdz [/(z)fg(z) 
Jdz [/(z)f • 



(8) 



This formulation has the virtue of making the q-expectation value of a constant equal to that constant. It 
has been found by Abe |l7j to resolve problems with the fmiteness of certain physical observables for the 
general-g ideal gas. It has also been found by Anteneodo |l8| to yield the same differential equation, albeit 
with renormalized coefficients, for the radial profile of the pure-electron plasma that was found in earlier 
studies [|]|7]] . In this work, however, we shall adhere to the original version of the second axiom. Most of 
the convergence problems encountered in the derivation of the equilibrium properties of the general-q ideal 
gas (lj] do not appear in our derivation of the hydrodynamic properties of that system. So, although the 
use of normalized g-expectation values would be an interesting modification to the current study, we leave 
it for future work. 



B. Global Hydrodynamic Equilibria 



We next consider phase space coordinates z = (x, v) , where x denotes position and v denotes velocity. 
Because our particles undergo only point collisions, we consider equilibria that conserve mass, momentum 
and kinetic energy. (We do not have to worry about the potential energy.) That is, we fix the g-expectation 
values of the n — 3 quantities 



7 1 (z) = m 
7 2 (z) = mv 
7 3 (z) = mv 2 /2, 



(9) 
(10) 
(11) 



namely 



/ dz [/(z)] 9 7 2 (z 



= / dz [f{z)Y 



mv 

mv 2 /2 



(12) 



The equilibrium distribution function is then 



/ (cq) (z) = 



1 + (q - 1)— Ai + A 2 • v + A 3 — 



(13) 



This may be written in the more concise form 



/ (eq) (z) =A 



l + (9-l) 



2k B T 



v — u 



(14) 



where the three quantities 



A=iq 



(15) 



A2 

"A 3 



(16) 
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T = 



and 

1 

^3 

are determined by the three requirements 
Af, 



in 



l + (,-l)-(A 



_ Af_\ 
2X 3 J 



(17) 



VA q / rfv 



1 + ^- 1 )€t' v - u ' 2 



mv 

mw 2 /2 



(18) 



and where V = J dx is the total spatial volume. The integrals can be done by making the substitution 

(v-u) 



w 



' m 
2k 



B 



9-1 



T 



so that 



M„ 



( 



mVA q 



2k B 
m 



T 



D/2 



J dw ( 



1 + aw 1 



\ 



u + 



2k B 


T 


V m 


9-1 





T 


u • w - 


k B 


T 




9-1 


m 


9-1 



(19) 

where we have defined a = sgn (T / (q — 1)) G {— 1, +1}- The region of integration in w is always spherically 
symmetric, so terms with odd integrands vanish, and we may use the Z?-dimensional spherically symmetric 
volume element 

27^ 



dw = 



r(f) 



) D - X dw 



(20) 



on the rest, obtaining 



M„ 



f 2nk B 



in 



2mVA q 
mVA q (2txI v . 



T 



m 



9-1 



T 



D/2 



dw w - 1 (l + aw 2 ) 



1 
u 



9-1 



D/2 



J dx x D ' 2 - 1 {l + axY 



u I k_B_ 

2 ^ m 



1 
U 



9-1 



9-1 



W 
\ 



(21) 



where x = w 2 . 

At this point, we need to specify the limits of integration. If a = +1, the integrand is defined for 
< x < oo. If a = — 1, then the integrand is defined only for < x < 1. In the latter case, states with x > 1 
or 



v-u > 



! 2k B 
m 



T 



9-1 



(22) 



are thermally forbidden. Thus, we have 



M„ 



mVA q f2nk B 



r(#) 



m 



T 



9-1 



D/2 



/ °° dxx D / 2 - 1 (l + xy^ 
dxx D ' 2 - 1 {l-x)-^ 



1 

u 

u 2 I kt 
2 ^ m 

1 

u 

u 2 I k B 
2 ' m 



9-1 



for a = +1 



9-1 



for (j 



(23) 
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In both cases, the integration results in a beta function, 
functions [fL3fl. The results are 



and this can be expressed in terms of gamma 



mVA q 



2nk B T 



D/2 



\i-i\ D/3 r{-^) 



/ 



u 

Dk B T/(2m) 



V 2 g-(?-i)(^+i) 
1 



u 



|l-<?r /2 r(T^ + -T) „ 2 , Dfc B T/(2m) 



for ct = +1 



for ct = — 1 



(24) 



2 1 l+(l-g)f 

In either case, we see that the momentum per unit mass, or hydrodynamic velocity, is given by 
The energy per unit mass is then 

= ^ q — Uq 
£q ~ M q ~ 2 +Lq ' 

where we have defined the thermal or internal energy per unit mass, 

-l 



Dk B T 
2m 



1 + (1 - 0) 



D 



(25) 



(26) 



(27) 



which, remarkably, is independent of a. 

At this point, we arrive at an interesting ambiguity We are tempted to interpret the variable that has been 
suggestively called T as the temperature of our physical system. Certainly, it reduces to the usual temperature 
for q = 1. If we make this interpretation, we might conclude that the familiar energy equipartition theorem 
is invalid for q ^ 1. That is, for general q the thermal energy per unit mass is apparently no longer directly 
proportional to the dimension D. Figs. (|l|, |J) plot the ratio i q j (DkBT/ (2m)) versus (q, D) for various values 
of (D, q), respectively. This ratio is equal to unity when equipartition holds. We see that it is less than unity 
for q < 1, and greater than unity for q > 1. We may even be tempted to attach physical significance to this 
result by noting that the equipartition theorem presumes ergodicity, and systems with q ^ 1 are generally 
not ergodic. This reasoning is faulty, however, because it depends upon our rather arbitrary definition of the 
temperature in terms of the Lagrange multipliers in Eq. (p^). If we had instead defined the temperature as 

T' = . T , „ , (28) 



1 + (!-«) 



which also reduces to the usual definition as q — > 1, we would have found that i q = DkgT' /(2m) and 
concluded that equipartition is valid for all q. The lesson here is that the definition of the temperature is 
rather arbitrary in thermodynamics, and we should not suppose that we can measure its absolute value in 
an experiment. So, even though we have theoretical reason to believe that t q {T) is g-dependent, we can not 
exploit this dependence to measure q experimentally To do that, it will be necessary to find two physical 
observables whose relationship to each other is g-dependent. We shall return to this point numerous times 
below. 

We also note from Eq. ( |2^ ) that there is a hard upper bound on q for this type of system. In order for the 
proportionality constant between i q and T not to become negative, it must be that 1 + (1 — q)D/2 > 0, or 

<?<! + §■ (29) 

This conclusion is independent of diffeomorphic transformations in our definition of T and therefore may 
be expected to have some general validity []. A similar inequality was derived in Q . We shall show below 



x Of course, it could well be that we want to generalize our notion of temperature to include negative values, in 
which case this upper bound no longer applies. 
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Energy Equipartition 
Versus q, Various D 




0.25 0.5 0.75 1 1.25 1.5 1.75 2 

q 

FIG. 1. A plot of the ratio t 9 /(fcsT/(2m)) versus q for D ranging from 1 to 4. The plot color becomes lighter as 
D increases. 



Energy Equipartition 
Versus D, Various q 




0.5 1 1.5 2 2.5 3 3.5 4 

D 

FIG. 2. A plot of the ratio t,/(iaT/(2m)) versus D for q ranging from 0.2 to 2.0 in increments of 0.2. The plot 
color becomes lighter as q increases. 
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Coefficient c Versus q, Various D 




FIG. 3. A plot of the coefficients c 9i r> versus q for D ranging from I to 4. The plot color becomes lighter as D 
increases. 



that hydrodynamic considerations - in particular, the positivity and finiteness of the thermal conductivity 
- impose an even more stringent upper bound on q. 



Finally, the normalization constant A can be expressed in terms of the mass density p q = -rr as follows 

D/2 I 

Pq 



A« = 



2nk B T 



for a = +1 



\i-Q\ D/2 r(r^+§) , 1 
r(X) «** = -l 



(30) 



It follows that the global hydrodynamic equilibrium (GHE) distribution function can be written in the form 

D/2~\ 



/M(z) = 
where we have defined 



Cq,DPq 

m 



2nk B T 



C„n 



k-H J/2 r(^ T ) 



k B T 



for a 



m ,2 
— v — u„ 

2 1 9 



(31) 



(32) 



for a = — 1, 



and with the proviso that /( cq )(z) = if the argument raised to the — 1/(<Z — 1) power in Eq. ( fjjj ) is 
negative. This is the generalization of the Maxwell-Boltzmann distribution function for the Generalized 
Thcrmostatistics. The coefficients c 9i d have a particularly simple form for even dimension D, 



D/2 

c q ,D = J] [i - (l - l)q] ; 



(33) 



in particular, we note that c 9i 2 = 1 and c q A — 2 — q. More generally, these coefficients arc plotted against q 
for various values of D in Fig. ||. Various distributions for D — 1 are illustrated in Figs. |J and S. 
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Generalized Distribution Function 
Fixed u and T, Varying q 



0.4 
0.3 

SZ-, 

I 0.2 
0.1 




FIG. 4. A plot of [/ (cq) (z)] q versus v for D = 1, p„ = 1, u, = 1/2, T = 1, and g ranging from 0.5 to 1.5 in 
increments of 0.1. The plot color becomes lighter as q increases. 



Generalized Distribution Function 
Fixed u and q, Varying T 
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0.3 




0.2 




0.1 









FIG. 5. A plot of the equilibrium distribution function versus v for D — 1, p q = 1, u q = 
ranging from 0.25 to 1.5 in increments of 0.25. The plot color becomes lighter as T increases. 



1/2, q = 1/2 and T 
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III. GENERALIZED KINETICS 

A. Local Hydrodynamic Equilibria 

Before introducing our generalization of the BGK collision operator, we first introduce the concept of 
local hydrodynamic equilibrium (LHE). The GHE distribution function derived in the preceeding section is 
a function of velocity v alone, and independent of position x. This is a direct consequence of the fact that 
the only conserved quantities present - mass, momentum and kinetic energy - are functions of v and not 
of x. The GHE distribution does depend parametrically on p q , u q and e q , so we may write it in the more 
precise format f^ eq \v; p q , u q: e q ). The LHE distribution function f^(z) is then defined as having the same 
functional form as the GHE, but with parameters p q , u q and e q weakly dependent on spatial position x. 
That is, 

/( )(z) = /(^(v ;/9g (x),u 9 (x),e g (x)). (34) 

It should be noted that the LHE, unlike the GHE, is not expected to be a stationary state of the system. 
If the system is initialized in a LHE distribution, its time evolution will generally take it away from this 
simple functional form. That is, the full solution to the system's kinetic equation should be expected to be 
the LHE plus a correction. To the extent that the spatial gradients are weak, this correction will be small 
and can be treated perturbatively. This is the basis of the Chapman-Enskog expansion. 



B. The Generalized BGK Equation and //-Theorem 



Armed with the concept of LHE, we propose the following generalization of the Boltzmann equation with 
the BGK collision operator: 



(35) 



where /^°^ (z) is the LHE distribution function that is defined as having the same moments as the exact 
distribution function /(z), 




dv [/(z)]< 



mv 

mv 2 /2 



dv 



/ (0) (z) 



mv 

mv 2 /2 



(36) 



and where r > is the collisional relaxation time. In spite of the deceptively simple appearance of Eq. (|3J[ 
note that it is nonlinear even when q = 1. This is because of Eq. (^) which requires that the moments 
of /(°)(z), which determine the parameters that appear nonlinearly in its functional form, are themselves 
functionals of /(z). To see this more clearly, note that Eqs. (34) and (|36|) can be used to rewrite Eq. ( ^B|) in 
the explicitly nonlinear but woefully cumbersome form 



d_ 

dt 



v-V [/(*)]« 



[/(z)P 



(eq) 



v; J dv [f{z)fm,J dv [f(z)]«mv, J dv [f(z)f 



(37) 



Hence, when q ^ 1, this departs from the usual BGK equation in two important respects: First, and most 
obviously, it is constructed so that it relaxes to the explicitly g-dependent LHE /(°)(z), given by Eqs. (|5|) and 
(pl|), rather than to the usual Maxwell-Boltzmann LHE. Second, the nonlinearity in the collision operator 
arises via the functional form of which also explicitly depends on q. It is important to keep these 

two mechanisms clearly in mind because a superficial glance at Eq. ([35]) might suggest that the substitution 
F(z) = Jf(z)] q will render its dynamics identical to those of the usual BGK equation, and this is not so. 

Eq. ( |35| ) is a sensible generalization of the BGK equation for three reasons. First, it obviously reduces to 
the usual BGK equation when q — 1. Second, if we multiply it by the 7*(z) and integrate over v, we find 
local conservation equations 
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mv 

TOW 
TOVW 2 /2 



(38) 



for each of the conserved quantities. Third, we can show that Eq. (|35|) obeys an ff-theorem - at least for 
q > - as follows: 



d 
If 



-Mfl 



dz 

~ T 

ks 
~ T 

ks 
~ T 



VWJ st 

_S 

St 



rfz 



1 



- i 

- i 



[/(*)]« 



rfz 



(g-l)r 

-g g [/ (0) ]-gg[/] 
r 

^[/ (0) ]-^[/] 

T 

5,[/W]-5 fl [/] 



1 



vim]*- 1 

k B 
(ff-l)r 
k B 



dz 



dzf 



(g-l)r 



+ V-VJ [/(*)]* 
{[/(*)]« -[/(°> (*)]«} 
/(^)-/ (0) (z) 
[<?(z) - 1] 



(0) 



[/(z)]g-[/(°>(z)]« 

[g(z)p-r 



(39) 



where g(z) = /(z)//(°)(z) > 0, the domain has been assumed to be without boundary, and where we have 
defined the function 



$,(*) = 



1 



(x-1) 



1 



9-1 



(40) 



The first term of Eq. ( p9| ) is nonnegative since the entropy is a maximum for the distribution function 
/(°)(z) and q > by construction. It is not difficult to show that the function <I> 9 is nonnegative for positive 
argument (see Fig. g), so the second term is also nonnegative. It follows that dS q [f]/dt > for q > 0, with 
equality holding only at equilibrium. 

For q < 0, the generalized entropy S q [f] is minimized at its extremum, rather than maximized, as can 
be seen from examination of its second variation. In this situation, the first term on the right-hand side of 
Eq. ( |39| ) is negative, but the second term is still positive because <& q (x) > still holds for q < 0. Hence the 
above argument breaks down and a more subtle one is needed. We shall not concern ourselves with the case 
q < in this paper. 



IV. THE CHAPMAN-ENSKOG ANALYSIS 



A. The Asymptotic Expansion 



We now develop the distribution function in a perturb ation series where the zero-order approximation is 
the LHE, following the program described in Subsection III A. We rewrite our generalized BGK equation, 
Eq. (p^), in the form 



(l + TT>)F(z)=FW(z), 



(41) 



where we have defined 
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0.5 1 1.5 2 2.5 3 3.5 4 

x 



FIG. 6. A plot of the function $ q (x) versus x for q 6 {0.25, 0.50, 0.75, 1.00, 1.25, 1.50}. The plot color becomes 
lighter as q increases. 



and 



and the operator 



F<°>(z) = [/(°)(z) 



(42) 

(43) 
(44) 



(45) 



The formal solution to the above equation is 

F(z) = (l + rV) F(°)(z). (46) 

Following the standard Chapman-Enskog procedure [ p"5| ] , we expand the derivative in a series of successively 
more slowly varying time scales, 



where 



where e is a formal expansion parameter, and t ? - is the jth time scale. To second order in e, we find 



F(z) = 1 - erVi - e 2 r (p 2 - rvfj +■■■ F (0) (z) 



(47) 



(48) 



(49) 



Thus we have 



F(z) = £\"FW(z), 

n=0 



(50) 
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where 



F (2) (z) = -T(D 2 -rDnFl u, (z) 



p(o), 



(51) 
(52) 



The higher order terms F^- 1 ' and -F^ 2 - 1 must not change the definition of the conserved densities, 



0= / dv 



F(z) - F (0) (z) ( mvv ] = e ' ' / dv ^ ( " } ( 
mvv 2 /2 



z) | mvv 

2 , 



mvt) /2 



(53) 



To fully specify the determination of the f( n ), we further require that they satisfy Eq. (|5^) order by order; 
that is, 



= f dv F (ll) (z) [ mvv ] 
J \ mvv 2 /2 J 



(54) 



for n > 1. 



B. The First-Order Solution 

At first order, Eq. ( pi} ) can be written 



i^(z) = -r(^- + v .v) F<°>(a). 



(55) 



We take the hydrodynamic moments of both sides, noting that the left-hand side will vanish thanks to 
Eq. (||). We get 



0=^| /M, 




mvv 

mvv 2 1 2 



The first of these is immediately seen to be the usual equation expressing conservation of mass 

= ^ + V-( Pq u q ), 



(56) 



(57) 



which is thus seen to be g-invariant at first order. 

To evaluate the momentum equation at first order, we must integrate the dyad mvv times F<W. We have 



/ dvF^(z) mvv = dv F^ (z)m [(v — u q ) + u q ] [(v — u 9 ) + u q 



p q u q u q + 



2k B 
m 



T 



dw F' '(z)raww, 



(58) 
(59) 



since the odd cross terms integrate to zero. Turning our attention to the last integral over w, it is clear from 
parity that only the diagonal elements of this dyad will have nonvanishing integral, and from isotropy that 
they will all give the same result. Thus, we have 



— +i 

JdvF^\ Z )mvv = Pq u q u q + l(^\-L T \y 2 J 



dw mu> . 



(60) 
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To perform this last integral, the spherically symmetric volume element of Eq. (|2C|) is inadequate. Rather, 
we write w x — w cos 9, where < 9 < 7r is a polar angle, and we adopt the cylindrically symmetric volume 
clement 



dw 



2?r— 

' D-l ^ 
s 2 , 



r 



-w 



D-X 



dw d9. 



(61) 



Now both the w and 9 integrals yield beta- functions which must be considered for both cases a = ±1. The 
result, which does not depend on a, is then 



J dv F (0) (z)ravv = p q u q u q + P q l 



where we have defined the pressure 



p q k B T 



m[l +(!-<?)§] 



The momentum equation at first order is thus 

d 



(62) 



(63) 



(64) 



Eq. (|63|) appears to be a nonideal equation of state for the pressure but, as noted at the end of the previous 
subsection, we must avoid attaching physical significance to the temperature T. Rather, we should take care 
to relate observable quantities - in this case, the pressure P q and the internal energy i q . From Eqs. (|63| ) and 
(p7|), we have 



Pq — jjPq L qi 



(65) 



and this is the usual ideal gas equation of state with no correction, which we see is g-invariant after all. 

Finally, we consider the energy equation at first order. We must integrate the vector mv» /2 times F^l 
We have 



dv F (o) {z) H^l = J dv F (o )(z) ™ [(v _ Ug) + Ug] | (v _ Ug) + Uq 



p q e q u q + u q 



2k B 
m 



T 



mww 



(66) 
(67) 



since the odd cross terms integrate to zero. The quantity in square brackets in the second term is the same 
one we encountered in the definition of the pressure and is equal to P q l. The energy equation at first order 
is thus 



d 

° = dt 1 ~ ^ Pq£q) + V ' ^ Pq£q + Pq) 



(68) 



Eqs. ([57j), ( |64| ) and ( |68|) give us the rate of change of the hydrodynamic variables on time scale t\ . They and 
the equation of state, Eq. (pq), are all seen to be q-invariant. Thus, at first order, there is no experiment that 
we could perform that would distinguish a system with q =/= 1 from the more orthodox q = 1 case. To make 
this distinction, we must examine the dissipative terms that appear at second order in the Chapman-Enskog 
analysis, and we turn our attention to those now. 



C. The Second-Order Solution 



At second order, Eq. (153) can be written 
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FW(z) = -T 



dt 2 



— +v.V)F(°)(z)-r( — 



dt 



-+vv) F (0) (z) 



(69) 



Again, we take the hydro-dynamic moments of both sides, noting that the left-hand side will vanish thanks 
to Eq. (|54|). Drawing on our experience from the first-order calculation, we see that we get 

PqUq 
PqU q U q +P q l 
{PqSg + Pq)Uq 





tV 



g I Pq u q \ r I " l 

— p q u q u q + P q l +V - / dvvvF(°'(z) mv 
dtl \ (p q e q + P q ) u q I J \ m« 2 /< 







/ m 


: / dvvvfC'fz) 


mv 






V mv 2 /2 


m \ 






mv 






mv 2 /2 J 







(70) 



where we used the first order results in the second step. Note that the right-hand side is now manifestly 
a divergence. The quantity inside this divergence is the negative of a diffusive flux. The derivatives with 
respect to t\ can be eliminated using the first-order evolution equations. 

The right-hand side of the second-order mass conservation equation is seen to vanish due to the first-order 
momentum conservation equation. Thus, we find that the form of the mass conservation equation 



0=^ £ + V-(p,u ? ) 



(71) 



is g-invariant to second order. 

To evaluate the right-hand side of the momentum equation at second order, we must integrate the triad 
ravw times F^°\ The procedure is the same as that at first order, and no new moments need to be defined. 
The result in index notation is 



dv (z)mViVjV k = p q u qi u qj u qk + P q (s i:j u qk + S ik u gj + S jk u q ^ . 



(72) 



Inserting this into the second-order momentum equation, carrying out the derivatives with respect to t\ 
using the first-order equations, and simplifying, we get 



d r 

— (p q U q ) + V ■ {p q U q U q + P q l) = V • | 



/<--/ 



Vu, + (Vu,f +A 9 (V-u 9 )l}, 



where the superscript T denotes "transpose," and where we have defined the shear viscosity 

2t 



and the bulk viscosity 



Pq — jj Pq l qi 



4r 



(73) 



(74) 



(75) 



Here we have used the equation of state, Eq. (|6q ) , to express the viscosities in terms of the internal energy. 
Note that Eqs. ([74]) and ( |75| ) are g-invariant, as is the ratio 



_£ 

Pq 



2 
D 



(76) 



which is known to be approximately true for a variety of real gases [ fL9[ . Thus, measurement of the ratio of 
the viscosities is still insufficient to determine q experimentally. 
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Finally, wc turn our attention to the energy equation at second order. For this we must integrate the 
diad mvw /2 times F^l The by now familiar procedure again yields beta-function integrals which must 
be considered for both cases a = ±1. The result, which does not depend on tr, is then 



dv F'°'(z)mvv- 



(p q e q + 2P q ) u q u q 



u\ 4 



1 + (1 - q) 



l + + 1) 



1. (77) 



Inserting this into the second-order energy equation, carrying out the derivatives with respect to t\ using 
the first-order equations, and simplifying, we get 



d 

— (p q e q ) + V • [( Pq e q + Pg) u q ] 



v • <! u, • 



(VU, + (VU/) + Xg (V • Ug) 1 



where we have defined the thermal conductivity 



k q Vi q 



(l-g)-*V0v,; 
Pq 



(78) 



k = — 

and the anomalous transport coefficient 



l + + 1) 



T Pq L 



(79) 



(I 



4 

q = ^2 



D 
2 



l + (l-g)(f + 1) 



Tp q L q . 



(80) 



We note that positivity and finiteness of the thermal conductivity sets an even more stringent upper bound 
on q than Eq. (E9), namely 



q<l 



D 



(81) 



We emphasize that this inequality is not expected to hold for general systems; its derivation was specific to 
our assumption of an ideal gas. 

The most striking feature of the second-order energy conservation equation, Eq. (|7S|), is that it is not 
g-invariant. Its diffusive flux contains a term proportional to (1 — q) V {p q L q ) / p q with a new transport 
coefficient a q . The presence of this term can be detected by purely hydrodynamic experiments and this may 
be used to test whether or not q is equal to unity. Moreover, we see that the ratio 



i + (i-?) T 



gives us another means of determining q, as does the ratio 



(82) 



/'-</ 



2 
D 



1 + (!-«) 



l + (l-g)(f + 1) 



(83) 



These ratios are plotted in Figs. [?] and | 



V. DISCUSSION 

Eqs. (|7l|), ( |73| ) and ( |78| ) constitute the compressible Navier-Stokes equations for the generalized thermo- 
statistics. We see that the mass and momentum conservation equations are completely q- invariant, but the 
energy conservation equation is not. It contains a term that is proportional to (1 — q) V (p q t q ) / p q which van- 
ishes for Boltzmann-Gibbs statistics. It is a term not normally included in presentations of the Navier-Stokes 
equations. 
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10 



Ratio a /k 

q q 

Versus q, Various D 




FIG. 7. A plot of the ratio of the anomalous transport coefficient a q to the thermal conductivity k q for D ranging 
from 1 to 4. The plot color becomes lighter as D increases. 



Ratio k ii 

q q 



Versus q, Various D 




1.75 



FIG. 8. A plot of the ratio of the thermal conductivity k q to the shear viscosity \x q for D ranging from 1 to 4. The 
plot color becomes lighter as D increases. 
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In addition, we have shown that certain ratios of the transport coefficients may be q dependent. While the 
famous result X q /n q = —2/D appears to be q- invariant, the ratios k q j [i q and a q /k q arc decidedly less robust 
and may be used to infer a value of q. At this point, the reader may be tempted to pick up a handbook of 
material properties, search for materials with anomalous ratios of thermal conductivity to shear viscosity, 
and attribute those anomalies to a breakdown in Boltzmann-Gibbs statistics for those materials. The reader 
is implored to resist this temptation. There are very many potential pitfalls in the Boltzmann/Chapman- 
Enskog analysis that are far more likely to generate such anomalies than a breakdown in the foundations 
of thermostatistics. The presence of rotational and/or internal molecular degrees of freedom, high densities 
leading to three-body collisions, and violations of the Boltzmann molecular chaos assumption are but three 
examples of phenomena that might also cause such anomalies. If, however, one has other reasons to believe 
that a particular system violates Boltzmann-Gibbs thermostatistics, then examination of the ratios of the 
transport coefficients may provide further corroboration. One does have such reason to believe this for 
the stellar polytrope and pure-electron plasma examples mentioned earlier; unfortunately, those systems of 
particles are not ideal gases in any reasonable sense of the term so, even if one could measure the transport 
coefficients of such things, the quantitative results derived above should not be expected to apply. 

In summary, the presence of the anomalous term in the energy equation of a dilute gas would seem to be 
the easiest way to search for deviations of q from unity. I am currently unaware of any experimental work 
that would clearly establish either the presence or absence of this term. I leave it to my colleagues who have 
more familiarity with the experimental hydrodynamic literature to sort out this matter. 



VI. CONCLUSIONS 



We have shown that the Navier-Stokes equations for mass and momentum conservation are g-invariant 
to second order in the Chapman-Enskog expansion, but that the equation for energy conservation is not. 
We have derived the form of this anomalous term, using a generalized Chapman-Enskog analysis on a 
generalized BGK kinetic equation. In addition, we have found g-dependent anomalies in ratios of certain 
transport coefficients - in particular, in the ratio of the thermal conductivity to the shear viscosity. Finally, 
we have found that hydrodynamics gives rise to the upper bound q < 1 + 2/(D + 2) which is more stringent 
than that imposed by equilibrium considerations. 

This work could be substantially improved by using a more realistic kinetic equation, but we have argued 
that the principal results - the form of the hydrodynamic equations, and the ratios of the transport coefficients 
- ought to be robust in this regard. It is hoped that some of the insights obtained in this analysis will 
eventually be useful in constructing simple experimental tests for the presence of breakdowns in Boltzmann- 
Gibbs statistics. 
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